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Abstract 

Numerous codes are being developed to solve Shallow Water equa- 
tions. Because there are used in hydraulic and environmental studies, 
their capability to simulate properly flow dynamics is critical to guaran- 
tee infrastructure and human safety. While validating these codes is an 
important issue, code validations are currently restricted because analytic 
solutions to the Shallow Water equations are rare and have been published 
on an individual basis over a period of more than five decades. This article 
aims at making analytic solutions to the Shallow Water equations easily 
available to code developers and users. It compiles a significant number of 
analytic solutions to the Shallow Water equations that are currently scat- 
tered through the literature of various scientific disciplines. The analytic 
solutions are described in a unified formalism to make a consistent set 
of test cases. These analytic solutions encompass a wide variety of flow 
conditions (supercritical, subcritical, shock, etc.), in 1 or 2 space dimen- 
sions, with or without rain and soil friction, for transitory flow or steady 
state. The corresponding source codes are made available to the commu- 
nity (http://www.univ-orleans.fr/mapmo/soft/SWASHES), so that users 
of Shallow Water-based models can easily flnd an adaptable benchmark 
library to validate their numerical methods. 
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1 Introduction 



Shallow- Water equations have been proposed by Adhcmar Barrc de Saint- Venant 
in 1871 to model flows in a channel [4]. Nowadays, they are widely used to model 
flows in various contexts, such as: overland flow [22, 47], rivers [25, 9], flood- 
ing [10, 17], dam breaks [1, 51], nearshore [6, 37], tsunami [23, 33, 41]. These 
equations consist in a nonlinear system of partial differential equations (PDE- 
s), more precisely conservation laws describing the evolution of the height and 
mean velocity of the fluid. 

In real situations (realistic geometry, sharp spatial or temporal variations 
of the parameters in the model, etc.), there is no hope to solve explicitly this 
system of PDE-s, i.e. to produce analytic formulae for the solutions. It is there- 
fore necessary to develop specific numerical methods to compute approximate 
solutions of such PDE-s, see for instance [49, 34, 7]. Implementation of any of 
such methods raises the question of the validation of the code. 

Validation is an essential step to check if a model (that is the equations, the 
numerical methods and their implementation) suitably describes the considered 
phenomena. There exists at least three complementary types of numerical tests 
to ensure a numerical code is relevant for the considered systems of equations. 
First, one can produce convergence or stability results {e.g. by refining the 
mesh). This validates only the numerical method and its implementation. Sec- 
ond, approximate solutions can be matched to analytic solutions available for 
some simplified or specific cases. Finally, numerical results can be compared 
with experimental data, produced indoor or outdoor. This step should be done 
after the previous two; it is the most difficult one and must be validated by a 
specialist of the domain. This paper focuses on the second approach. 

Analytic solutions seem underused in the validation of numerical codes, pos- 
sibly for the following reasons. First, each analytic solution has a limited scope 
in terms of fiow conditions. Second, they are currently scattered through the 
literature and, thus, are difficult to find. However, there exists a significant 
number of published analytic solutions that encompasses a wide range of flow 
conditions. Hence, this gives a large potential to analytic solutions for validata- 
tion of numerical codes. 

This work aims at overcoming these issues, on the one hand by gathering 
a significant set of analytic solutions, on the other hand by providing the cor- 
responding source codes. The present paper describes the analytic solutions 
together with some comments about their interest and use. The source codes 
are made freely available to the community through the SWASHES (Shallow- 
Water Analytic Solutions for Hydraulic and Environmental Studies) software. 
The SWASHES library does not pretend to list all available analytic solutions. 
On the contrary, it is open for extension and we take here the opportunity to ask 
users to contribute to the project by sending other analytic solutions together 
with the dedicated code. 

The paper is organized as follows: in Section 2, we briefiy present the nota- 
tions we use and the main properties of Shallow- Water equations. In the next 
two sections, we briefiy outline each analytic solution. A short description of the 
SWASHES software can be find in Section 5. The final section is an illustration 
using the results of the Shallow- Water code developed by our team for a subset 
of analytic solutions. 
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Figure 1: Notations for 2D Shallow- Water equations 



2 Equations, notations and properties 

First we describe the rather general settings of viscous Shallow- Water equations 
in two space dimensions, with topography, rain, infiltration and soil friction. 
In the second paragraph, we give the simplified system arising in one space 
dimension and recall several classical properties of the equations. 

2.1 General settings 

The unknowns of the equations are the water height {h{t, x, y) [L]) and u{t, x, y), 
v{t,x,y) the horizontal components of the vertically averaged velocity [L/T] 
(Figure 1). 

The equations take the following form of balance laws, where g ~ 9.81 m/s^ 
is the gravity constant: 

dth + dx (hu) + dy (hv) = R-I 

dt {hu) + dx (^hu'^ + + dy (huv) = gh{Sox - SfJ + fj,Sdx ^^-^ 

dt (hv) + dx (huv) + dy (^hv^ + = gHSoy - Sjy) + fiSdy 

The first equation is actually a mass balance. The fiuid density can be replaced 
by the height because of incompressibility. The other two equations are mo- 
mentum balances, and involve forces such as gravity and friction. We give now 
a short description of all the terms involved, recalling the physical dimensions. 

• z is the topography [L], since we consider no erosion here, it is a fixed 
function of space, z{x, y), and we classically denote by Sqx (resp. S^y) the 
opposite of the slope in the x (resp. y) direction, = —dxz{x,y) (resp. 
Soy = -dyz{x,y))] 

• i? > is the rain intensity [L/T], it is a given function R{t,x,y) > 0. In 
this paper, it is considered uniform in space; 



I is the infiltration rate [L/T], mentioned for the sake of completeness. It 
is given by another model (such as Green-Ampt, Richards, etc.) and is 
not taken into account in the following; 



S f ~ \ S f ^fy) ^^"^ friction force. The friction law Sf may take several 
forms, depending on both soil and flow properties. In the formula below, 
U is the velocity vector U = {u,v) with \U\ = \/v? + v'^ and Q is the 
discharge Q = (hu,hv). In hydrological models, two families of friction 
laws are encountered, based on empirical considerations. On the one hand, 
we have the family of Manning-Strickler's friction laws 

_r -r ^1^1 

Cf = 71^, where n is the Manning's coefficient [L'^/^^T]. 

On the other hand, the laws of Darcy-Weisbach's and Chezy's family 

writes 

_r -r ^1^1 

With C/ = f/{8g), f dimensionless coefficient, (resp. Cj = C 
[L^/^/T]) we get the Darcy-Weisbach's (resp. Chezy's) friction law. Notice 
that the friction may depend on the space variable, especially on large 
parcels. In the sequel this will not be the case. 

finally, ^Sd ~ {^piSdx-, ^'Sdy) is the viscous term with /i > the viscosity 
of the fluid [LVT]. 



2.2 Properties 

In this section, we recall several properties of the Shallow- Water model that are 
useful in the flow description. For the sake of simplicity, we place ourselves in the 
one-dimensional case, extensions to the general setting being straightforward. 
The two-dimensional equations (1) rewrite 



dth + dx{hu) = R 



dtihu) + dx yiu^ + — j = gh{Sox - Sf) + fid^ (hd^u) 

The left-hand side of this system is the transport operator, corresponding to 
the flow of an ideal fluid in a flat channel, without friction, rain or infiltration. 
This is actually the model introduced by Saint- Venant in [4], and it contains 
several important properties of the flow. In order to emphasize these properties, 
we flrst rewrite the one-dimensional equations using vectors form: 

dtW + dxF{W) = 0, where W ^ j^^^ F{W) = ^ ^^2'^ 3^ j ' (3) 

with F{W) the flux of the equation. The transport is more clearly evidenced 
in the following nonconservative form, where yl(VF) — F'{W) is the matrix of 
transport coefficients: 

dtW + A{W)dxW = 0, A{W) = F'{W) ^l^^2lg^ 2I ) • 
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More precisely, when h > 0, the matrix turns out to be diagonalizable, 

with eigenvalues 

Xi{W) = u - <u+ yfgh = A2(W). 

This important property is called strict hypcrbolicity (sec for instance [24] and 
references therein for more complete information). The eigenvalues are indeed 
velocities, namely the ones of surface waves on the fluid, which are basic char- 
acteristics of the flow. Notice here that the eigenvalues coincide if /i = m, that 
is for dry zones. In that case, the system is no longer hyperbolic, and this in- 
duces difficulties at both theoretical and numerical levels. Designing numerical 
schemes that preserve positivity for h is very important in this context. 

From these formula we recover a useful classification of flows, based on the 
relative values of the velocities of the fiuid, u, and of the waves, ^/gh. Indeed 
if |u| < \/gh the characteristic velocities have opposite signs, and information 
propagate upward as well as downward the flow, which is then said subcritical 
or fluvial. On the other hand, when |m| > ^/gh, the flow is supercritical, or 
torrential, all the information go downwards. A transcritical regime exists when 
some parts of a flow are subcritical, other supercritical. 

Since we have two unknowns h and u (or equivalently h and q = hu), a 
subcritical flow is therefore determined by one upstream and one downstream 
value, whereas a supercritical flow is completely determined by the two upstream 
values. Thus for numerical simulations, we have to impose one variable for 
subcritical inflow/outflow. We impose both variables for supercritical inflow 
and for supercritical outflow, free boundary conditions are considered (sec for 
example [8]). 

In this context, two quantities are useful. The flrst one is a dimcnsionlcss 
parameter called the Froude number 

Fr^^. (5) 
y/gh 



It is the analogue of the Mach number in gas dynamics, and the flow is subcritical 
(resp. supercritical) if Fr < 1 (resp. Fr > 1). The other important quantity is 
the so-called critical height he which writes 



/ \ 2/3 



(6) 



for a given discharge q ~ hu. It is a very readable criterion for criticality: the 
flow is subcritical (resp. supercritical) ii h > he (resp. h < he). 

When additional terms are present, other properties have to be considered, 
for instance the occurrence of steady states (or equilibrium) solutions. These 
specific flows are defined and discussed in Section 3. 

3 Steady state solutions 

In this section, we focus on a family of steady state solutions, that is solutions 
that satisfy: 

dth — dtu — 0. 
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Replacing this relation in the one-dimensional Shallow- Water equations (2), the 
mass equation gives dx(hu) = R or hu = q = Rx + qo, where go = liti^ = 0)- 
Similarly the momentum equation writes 

dx + = -ghdxZ - ghSf{h,q) + pidx (j^'^^j^ 

Thus for /i 7^ 0, we have the following system 

q^ Rx + qo, 

dx.^l^[^-gh)dxk-S,ik,q) + -^^Ox{hdxl). 

System (7) is the key point of the following series of analytic solutions. For 
these solutions, the strategy consists in choosing either a topography and get- 
ting the associated water height or a water height and deducing the associated 
topography. 

Since [5], it is well known that the source term treatment is a crucial point 
in preserving steady states. With the following steady states solutions, one can 
check if the steady state at rest and dynamic steady states are satisfied by the 
considered schemes using various fiow conditions (fluvial, torrential, transcriti- 
cal, with shock, etc.). Moreover, the variety of inflow and outflow configurations 
(flat bottom/ varying topography, with/without friction, etc.) gives a validation 
of boundary conditions treatment. One must note that, as different source terms 
(topography, friction, rain and diffusion) arc taken into account, these solutions 
can also validate the source terms treatment. 

The last remark deals with initial conditions: if initial conditions are taken equal 
to the solution at the steady state, one can only conclude on the ability of the 
numerical scheme to preserve steady states. In order to prove the capacity to 
catch these states, initial conditions should be different from the steady state. 
This is the reason why the initial conditions, as well as the boundary conditions, 
arc described in each case. 

Table 1 lists all steady-state solutions available in SWASHES and outlines 
their main features. 



3.1 Bumps 

Here we present a series of steady state cases proposed in [27, p. 14-17] based on 
an idea introduced in [32], with a flat topography at the boundaries, no rain, 
no friction and no diffusion (i? = m/s, Sf — and ^ ~ m^/s). Thus system 
(7) reduces to 

' q = qa, ^ 

In the case of a regular solution, we get the Bernoulli relation 

^-HM-) + ^(x) = C.t (8) 
which gives us the link between the topography and the water height. 
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steady- state solutions 


Flow criticality 


Friction 




Type 


Description 


§ 


Reference 


Sub. 


Sup. 


Sub . ^ Sup. 


Jump 


Man. 


D.-W. 


Other 


Null 


Comments 


Bumps 


Lake at rest with, immersed bump 


3.1.1 


[16] 
















X 


Hydrostatic equilibria 


Lake at rest with emerged bump 


3.1.2 


[16] 
















X 


Hydrostatic equilibria and 
wet-dry transition 


Subcritical flow 


3.1.3 




27 




X 














X 


Initially steady state at rest 


Transcritieal flow without shock 


3.1.4 




27 








X 










X 


Initially steady state at rest 


Transcritical flow with shock 


3.1.5 




27 










X 








X 


Initially steady state at rest 


Flumes 
(MacDonald's 
based) 


with subcritical flow 


3.2.1 


[53] 


X 








X 


X 






Initially dry channel. 1000 m long 


with supercritical flow 


3.2.1 


[16] 




X 






X 


X 






Initially dry channel. 1000 m long 


with sub- to super-critical flow 


3.2.1 


[53] 






X 




X 


X 






Initially dry channel. 1000 m long 


with super- to sub-critical flow 


3.2.1 


[53] 








X 


X 


X 






Initially dry channel. 1000 m long 


Short channel 
with smooth transition and shock 


3.2.2 


[53] 






X 


X 


X 








At t=0. lake downstream. 
100 m long 


bikort channel 
with supercritical flow 


3.2.2 


[16] 




X 






X 








Initially dry. 100 m long 


Short channel 
with sub- to super-critical flow 


3.2.2 


[53] 






X 




X 








At t=0, lake downstream. 
100 m long 


Very long, undulating and periodic channel 
with subcritical flow 


3.2.3 


[53] 


X 








X 








At t=0, lake downstream. 
5000 m long 


Rain on a long channel 
with subcritical flow 


3.3.1 


[53] 


X 








X 


X 






Initially dry. 1000 m long 


Rain on a long channel 
with supercritical flow 


3.3.2 


[53] 




X 






X 


X 






Initially dry. 1000 m long 


with subcritical flow and diffusion 


3.4.1 


[19] 


X 












X 




Initially dry. 1000 m long 


with supercritical flow and diffusion 


3.4.2 


[19] 




X 










X 




Initially dry. 1000 m long 


PsGudo-2D short channel 
with subcritical flow 


3.5.1 


[35] 


X 








X 








Initially dry. 200 m long 


PsGudo-2D short channel 
with supercritical flow 


3.5.2 


[35] 




X 






X 








Rectangular cross section. 
Initially wet. 200 m long 


Pscudo-2D short channel 
with smooth transition 


3.5.3 


[35] 






X 




X 








Rectangular cross section. 
Initially part!y-wet. 200 m long 


Pscudo-^D short channel 
with shock 


3.5.4 


[35] 








X 


X 








Rectangular cross section. 
Initially dry. 200 m long 


Pscudo-^D long channel 
with subcritical flow 


3.5.5 


[35] 


X 








X 








Initially dry. 400 m long 


Pscudo-2D long channel 
with smooth transition and shock 


3.5.6 


[35] 






X 


X 


X 








Isoscclc trapezoidal cross section. 
Initially dry. 400 m long 



Slope is always variable. Sub.: Subcritical; Sup.: Supercritical; Man.: Manning: D.-W.: Darcy-Weisbach 



Table 1: Analytic solutions for shallow flow equations and their main features 
— Steady-state cases 
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Initial conditions satisfy the hydrostatic equihbrium 

h + z ^ Cst and q^O m^/s. (9) 

These solutions test the preservation of steady states and the boundary condi- 
tions treatment. 

In the following cases, we choose a domain of length L = 25 m with a 
topography given by: 

, , J 0.2 - 0.05(a; - 10)^ if 8 m < .t < 12 m, 
^^'^> " \ else. 

3.1.1 Lake at rest with an immersed bump 

In the case of a lake at rest with an immersed bump, the water height is such that 
the topography is totally immersed [16]. In such a configuration, starting from 
the steady state, the velocity must be null and the water surface should stay flat. 

In SWASHES we have the following initial conditions: 

h + z = 0.5 m and q ~ m^/s 
and the boundary conditions 

r h = 0.5m, 
1^ = m^/s. 

3.1.2 Lake at rest with an emerged bump 

The case of a lake at rest with an emerged bump is the same as in the previ- 
ous section except that the water height is smaller in order to have emergence 
of some parts of the topography [16]. Here again, we initialize the solution at 
steady state and the solution is null velocity and flat water surface. 

In SWASHES we consider the following initial conditions: 
h + z ^ max(0.1, z) m and q — m^/s 
and the boundary conditions 

r ft. = 0.1m, 
[ g = m^/s. 

3.1.3 Subcritical flow 

After testing the two steady states at rest, the user can increase the difficulty 
with dynamical steady states. In the case of a subcritical flow, using (8), the 
water height is given by the resolution of 

h{xf + (^z{x) -^^-hi}j h{xf + |i = 0, Vx e [0,i], 

where = h{x ~ L). 
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Thanks to the BernouUi relation (8), we can notice that the water height 
is constant when the topography is constant, decreases (respectively increases) 
when the bed slope increases (resp. decreases). The water height reaches its 
minimum at the top of the bump [27]. 

In SWASHES, the initial conditions are 

h + z = 2 m and q = /s 

and the boundary conditions are chosen as 

upstream: q = 4.42 m^/s, 
downstream: h ~ 2 m. 

3.1.4 Transcritical flow without shock 

In this part, we consider the case of a transcritical flow, without shock [27]. 
Again thanks to (8), we can express the water height as the solution of 

h{xf + (zix) ~^-h,~ zm) h{xf + ^ = 0, Vx e [0, L], 

where zm = ma-x^Gfo.L] z and he is the corresponding water height. 

The flow is fluvial upstream and becomes torrential at the top of the bump. 
Initial conditions can be taken equal to 

h + z = 0.66 ni and q — m^/s 

and for the boundary conditions 

upstream: q = 1.53 m^/s, 

downstream: h = 0.66 m while the flow is subcritical. 

3.1.5 Transcritical flow with shock 

If there is a shock in the solution [27], using (8), the water height is given by 
the resolution of 

h{XshockY + \ z{Xshock) - - he - Z},[ \ h{Xshoek)'^ ~^ 



h{XshockY + ( z{Xshock) - - hi}j h{Xshockf + 1^ ^ ^^^^ 



In these equalities, zm = max^g^Q^^j z, he is the corresponding water height, 
/iL = h{x = L) and hi, ft,2 are the water heights upstream and downstream 
respectively. The shock is located thanks to the third relation in system (10), 
which is a Rankine-Hugoniot's relation. 

As for the previous case, the flow becomes supercritical at the top of the 
bump but it becomes again fluvial after a hydraulic jump. 
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One can choose for initial conditions 



h + z = 0.33 m and q = m^/s 



and the foUowing boundary conditions 



{ 



upstream: q ~ 0.18 m^/s, 
downstream: h ~ 0.33 m. 



We can find a generahsation of this case with a friction term in Hervouet's 
work [31]. 

3.2 Mac Donald's type ID solutions 

Following the lines of [35, 36], we give here some steady state solutions of sys- 
tem (2) with varying topography and friction term (from [53, 16]). Rain and 
diffusion are not considered {R ~ m/s, ^ — m^/s), so the steady states 
system (7) reduces to 



From this relation, one can make as many solutions as required. In this section, 
we present sonic of them that are obtained for specific values of the length L 
of the domain, and for fixed parameters (such as the friction law and its coef- 
ficient). The water height profile and the discharge are given, and we compute 
the corresponding topographies solving equation (11). We have to mention that 
there exists another approach, classical in hydraulics. It consists in considering 
a given topography and a discharge. From these, the steady-state water height 
is deduced thanks to equation (11) and to the classification of water surface pro- 
files (see among others [14] and [30]). Solutions obtained using this approach 
may be found in [56] for example. Finally, note that some simple choices for the 
free surface may lead to exact analytic solutions. 

The solutions given in this section are more intricate than the ones of the 
previous section, as the topography can vary near the boundary. Consequently 
they give a better validation of the boundary conditions. If 5/ 7^ (we have 
friction at the bottom), the following solutions can prove if the friction terms 
are coded in order to satisfy the steady states. 

Remark 1 All these solutions are given by the numerical resolution of an equa- 
tion. So, the space step should be small enough to have a sufficiently precise so- 
lution. It means that the space step used to get these solutions should be smaller 
than the space step of the code to be validated. 

3.2.1 Long channel: 1000 m 

Subcritical case We consider a 1000 m long channel with a discharge of 
q — 2 m^/s [53]. The flow is constant at infiow and the water height is prescribed 
at outfiow, with the following values: 




(11) 



upstream: q = 2 m^/s, 
downstream: h = hexi 
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The channel is initially dry, i.e. initial conditions are 

= m and q = Q vc? /&. 
The water height is given by 

'...(x)^(l)"'(l + lo.p(- 16(^-1)")), ,12) 

We remind the reader that q = 2 ni^/s on the domain and that the topography 
is calculated iteratively thanks to (11). We can consider the two friction laws 
explained in the introduction, with the coefficients n = 0.033 mr^/^s for Man- 
ning's and / = 0.093 for Darcy-Weisbach's. 
Under such conditions, we get a subcritical steady flow. 



Supercritical case We still consider a 1000 m long channel, but with a con- 
stant discharge q = 2.5 m^/s on the whole domain [16]. The flow is supercritical 
both at inflow and at outflow, thus we consider the following boundary condi- 
tions: 

J upstream: g = 2.5 m^/s and h = hex{Q), 
\ downstream: free. 

The initial conditions arc a dry channel 

h — m and q = m^/s. 
With the water height given by 




and the friction coefficients equal to n = 0.04 m"^/^s for Manning's and / = 
0.065 for Darcy-Weisbach's friction law, the flow is supercritical. 



Subcritical-to-supercritical case The channel is 1000 m long and the dis- 
charge at equilibrium is q ~ 2 m^/s [53]. The flow is subcritical upstream and 
supercritical downstream, thus we consider the following boundary conditions: 

J upstream: q ^ 2 m^/s, 
1 downstream: free. 

As initial conditions, wc consider a dry channel 

h = m and q = m^/s. 

In this configuration, the water height is 



^e2:(^) — 




for m < x < 500 m, 
for 500 m < a; < 1000 m. 



with a friction coefficient n = 0.0218 m"^/'^s (resp. / = 0.042) for the Manning's 
(resp. the Darcy-Weisbach's) law. Thus we get a transcritical flow (from fluvial 
to torrential via a transonic point). 
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Supercritical-to-subcritical case As in the previous cases, the domain is 
1000 ni long and the discharge is q = 2 m^/s [53]. The boundary conditions are 
a torrential inflow and a fluvial outflow: 

J upstream; q ~ 2 m^/s and h — h^xiO), 
y downstream: h = hexiWOO). 

At time t — s, the channel is initially dry 

h — ni and q — m^/s. 



The water height is defined by the following discontinuous function 




for 500 m < a; < 1000 m, 
with ai = -0.348427, aa = 0.552264, as = -0.55558. 

The friction coefficients are n = 0.0218 m'^^^s for the Manning's law and / = 
0.0425 for the Darcy-Weisbach's law. The steady state solution is supercritical 
upstream and becomes subcritical through a hydraulic jump located at a; = 



500 m. 

3.2.2 Short channel: 100 m 

In this part, the friction law we consider is the Manning's law. Generalization 
to other classical friction laws is straightforward. 

Case with smooth transition and shock The length of the channel is 100 
m and the discharge at steady states is q ~ 2 m^/s [53]. The flow is fluvial both 
upstream and downstream, the boundary conditions are flxed as follows 

J upstream: q = 2 m^/s, 

1 downstream: h = hexi^OO)- 

To have a case including two kinds of flow (subcritical and supercritical) and 
two kinds of transition (transonic and shock), we consider a channel filled with 
water, i.e. 

h{x) = max(/ie:,(100) + z(lOO) - z{x),Q) and g = m^s. 
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The water height function has 



the following formula 



hex {^) 



1/3 



1/3 



4 X \ 
3 ~ 100 / 



9x 
1000 



X 

100 



for m < a; < 




66.67 m, 



ai 



X 2 

Too " 3 

2 

^ 3 



a3 



ai 



Vloo 

2 
3 



04 



Vloo 

, 200 

for — — « 66.67 m < a; < 100 m, 



with fli = 0.674202, 02 = 21.7112, 03 = 14.492 et 04 = 1.4305. 
In this case, the Manning's friction coefficient is n = 0.0328 m-i/^s, the inflow 
is subcritical, becomes supercritical via a sonic point, and, through a shock 
(located at x = 200/3 « 66.67 m), becomes subcritical again. 



Supercritical case The channel we consider is still 100 m long and the equi- 
librium discharge is g = 2 m^/s [16]. The flow is torrential at the bounds of the 
channel, thus the boundary conditions are 

J upstream: q = 2 m^/s and h = hex{0), 
1^ downstream: free. 

As initial conditions, we consider an empty channel which writes 

h = m and q = m^/s. 

The water height is given by 




and the friction coefficient is n = 0.03 m'^/'^s (for the Manning's law). The flow 
is entirely torrential. 



Subcritical-to-supercritical case A 100 m long channel has a discharge of 
q = 2 m^/s [53]. The flow is fluvial at inflow and torrential at outflow with 
following boundary conditions 

J upstream: q ^2 m^/s, 
1 downstream: free. 

As in the subcritical case, the initial condition is an empty channel with a puddle 
downstream 

h{x) = max(/ie:c(100) + ;z(100) - z{x),0) and q = m^/s, 
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and the water height is 



The Manning's friction coefficient for the channel is n = 0.0328 m-i/3s. We get 
a transcritical flow: subcritical upstream and supercritical downstream. 

3.2.3 Long channel: 5000 m, periodic and subcritical 

For this case, the channel is much longer than for the previous ones: 5000 m, 
but the discharge at equilibrium is still g = 2 m^/s [53]. Inflow and outflow are 
both subcritical. The boundary conditions are taken as: 

J upstream: q = 2 m^/s, 

1 downstream: ft, = /lea; (5000). 

We consider a dry channel with a little lake at rest downstream as initial con- 
ditions: 

h{x) = max(/ie;j;(5000) + z(5000) - z{x),0) and g = mVs. 
We take the water height at equilibrium as a periodic function in space, namely 

, , , 9 1 . / TTX \ 

hpTix) = — \ — sm 

^ ' 8 4 V500/ 

and the Manning's constant is n = 0.03 m'^/'^s. We get a subcritical flow. As 
the water height is periodic, the associated topography (solution of Equation 
(11)) is periodic as well: we get a periodic configuration closed to the ridges- 
and-furrows configuration. Thus this case is interesting for the validation of 
numerical methods for overland flow simulations on agricultural fields. 

3.3 Mac Donald's type ID solutions with rain 

In this section, we consider the Shallow Water system (2) with rain (but without 
viscosity: fj, = m^/s) at steady states [53]. The rain intensity is constant, equal 
to Rq. The rain is uniform on the domain [0,ij]. Under these conditions, if we 
denote by go the discharge value at infiow q{t,0) = qq, we have: 

q{x) =qo + xRo, for < x < L. (14) 

The solutions are the same as in section 3.2, except that the discharge is given 
by (14). But the rain term modifies the expression of the topography through 
a new rain term as written in (7). More precisely. Equation (11) is replaced by 

f 9' 1^ ^ 2gi?o 

These solutions allow the validation of the numerical treatment of the rain. 
Remark 1, mentioned in the previous section, applies to these solutions too. 
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3.3.1 Subcritical case for a long channel 

For a 1000 m long channel, wc consider a flow which is fluvial on the whole 
domain. Thus wc impose the following boundary conditions: 

J upstream: q = qo, 

\ downstream: = /iex (1000), 

with the initial conditions 

h = m and q = m^/s, 

where hex is the water height at steady state given by (12). 

In SWASHES, for the friction term, we can choose either Manning's law 
with n = 0.033 m'^/'^s or Darcy-Weisbach's law with / = 0.093, the discharge 
qo is fixed at 1 m^/s and the rain intensity is Rq ~ 0.001 m/s. 



3.3.2 Supercritical case for a long channel 

The channel length remains unchanged (1000 m), but, as the flow is supercriti- 
cal, the boundary conditions are 

upstream: q ^ qo and h = hexiO), 
downstream: free. 

At initial time, the channel is dry 

ft, = m and q = m^/s. 

At steady state, the formula for the water height is (13). 

From a numerical point of view, for this case, we recommend that the rain 
does not start at the initial time. The general form of the recommended rainfall 
event is 

m/s if t < tji, 
Ro else, 



with tn = 1500 s. Indeed, this allows to get two successive steady states: for 
the first one the discharge is constant in space qo and for the second one the 
discharge is (14) with the chosen height profile (13). 

In SWASHES, we have a friction coefficient n = 0.04 m'^/'^s for Manning's 
law, / = 0.065 for Darcy-Weisbach's law. Inflow discharge is go = 2.5 m^/s and 
Ro = 0.001 m/s. 



3.4 Mac Donald's type ID solutions with diffusion 

Following the lines of [35, 36], Delestre and Marche, in [19], proposed new 
analytic solutions with a diffusion term (with R = m/s). To our knowledge, 
these are the only analytic solutions available in the litterature with a diffusion 
source term. In [19], the authors considered system (2) with the source terms 
derived in [38], i.e.: 

1 ( ao(h) /, M I \ , 

' -u + ai(n)\u\u \ and /i = 4/i/j 



9 \ h 
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with ^ ^ 

Oio(h) = T~r c(i{h) = — 

1 + 1^ fi + M 

where /it, [T] (respectively /i/j [L^/T]) is the vertical (resp. the horizontal) eddy 
viscosity and ki [T/L] (resp. kt [1/L]) the laminar (resp. the turbulent) friction 
coefBcicnt. At steady states we recover (7) with R = m/s, or again 

q = 90, 

d.z = i- - gh^ d^h - Sf{h, q) + -^ {-qd^.h + lid.h)^) . ^^^^ 

As for the previous cases, the topography is evaluated thanks to the momentum 
equation of (15). 

These cases allow for the validation of the diffusion source term treatment. 
These solutions may be easily adapted to Manning's and Darcy-Wcisbach's fric- 
tion terms. Remark 1 applies to these solutions too. 

In [19], the effect of /i^. fiy, kt and ki is studied by using several values. In 
what follows, we present only two of these solutions: a subcritical flow and a 
supercritical flow. 



3.4.1 Subcritical case for a long channel 

A 1000 m long channel has a discharge oi q ^ 1.5 m^/s. The flow is fluvial at 
both channel boundaries, thus the boundary conditions arc: 

J upstream: q = 1.5 m^/s, 

[ downstream: /i = /lea, (1000). 

The channel is initially dry 

h = m and q = m^/s 

and the water height at steady state (hex) is the same as in section 3.2.1. 

In SWASHES, the parameters are: h = 0.01, ki = 0.001, ^j.^ = 0.01 and 
tih = 0.001. 



3.4.2 Supercritical case for a long channel 

We still consider a 1000 m long channel, with a constant discharge q = 2.5 m^/s 
on the whole domain. Inflow and outflow arc both torrential, thus we choose 
the following boundary conditions: 

r upstream: q = 2.5 m^/s and h = hex{0), 
\ downstream: free. 

We consider a dry channel as initial condition: 

h ^ m and q = m^/s. 

The water height hex at steady state is given by function (13). 

In SWASHES, we have: h = 0.005, h = 0.001, = 0.01 and fih = 0.1. 
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3.5 Mac Donald pseudo-2D solutions 



In this section, we give several analytic solutions for the pseudo-2D Shallow- 
Water system. This system can be considered as an intermediate between the 
one-dimensional and the two-dimensional models. More precisely, these equa- 
tions model a flow in a rectilinear three-dimensional channel with the quantities 
averaged not only on the vertical direction but also on the width of the channel. 
For the derivation, see for example [26]. Remark 1, mentioned for Mac Donald's 
type ID solutions, applies to these pseudo-2D solutions too. 

We consider six cases for non-prismatic channels introduced in [35]. These 
channels have a variable slope and their width is also variable in space. More 
precisely, each channel is determined through the definition of the bottom width 
B (as a function of the space variable x) and the slope of the boundary Z 
(Figure 2). The bed slope is an explicit function of the water height, detailed 
in the following. 




Figure 2: Notations for the Mac Donald pscudo-2D solutions 



Solution 


B (m) 


Z (m) 


i(m) 


hin (m) 


hout (m) 


Subcritical flow in a short 
domain 


Biix) 





200 




0.902921 


Supercritical flow in a short 
domain 


Bi{x) 





200 


0.503369 




Smooth transition in a short 
domain 


Bi{x) 





200 






Hydraulic jump in a short 
domain 


Bi{x) 





200 


0.7 


1.215485 


Subcritical flow in a long 
domain 


B2{x) 


2 


400 




0.904094 


Smooth transition followed 
by a hydraulic jump in a 
long domain 


B2{x) 


2 


400 




1.2 



Table 2: Main features of the cases of pseudo-2D channels 
The features of these cases are summarized in table 2. In this table, the functions 
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B(x)/2 and -B(\)/2 for B=B1 



B(xV2 and -B(x)/2 for B=B2 




Figure 3: Functions Bi and B2 for tlie sliape of the channel 



for the bed shape are: 



Bi{x) = 10 - 5exp ( -10 - i j j for m < x < L = 200 m 

X 1\\ I X 2 



2N 



B.i.) = 10-5 exp (^-50 - 3 j j - 5 exp (^-50 ^- - , 

for m < X < L = 400 m , 
(Figure 3) and hm (resp. hout) is the water heiglit at the inflow (resp. outflow). 

In each case, the Manning's friction coefficient is n = 0.03 ni"^/^s, the discharge 
q is taken equal to 20 m'^s"^, the slope of the topography is given by: 

q^{B{x)+2Zh{x)) ^ • 



S^ix) = 1 - 



V h{xY^/^iBix) + Zh(x))^ 



gh{xf{B{x) + Zh{x)YJ h{xy"/3 {B{x) + Zh{x))^°'^ 

q^B'jx) 

~ gh{xY(B{x) + Zh{x) f 

where h is the water height and the topography is defined as z{x) — / So{X) dX. 

J X 

Remark 2 When programming these formulae, we noted a few typos in [35], 
in the expression of So{x), of 4> in the hydraulic jump case ('§ 3.5.4) o.iT'd in the 
formulation of h in [0; 120] in the solution for the smooth transition followed by 
a hydraulic jump 3.5.6). 

Remark 3 We recall that the following analytic solutions are solutions of the 
pseudo-2D Shallow-Water system. This is the reason why, in this section, h^x 
does not depend on y. 



3.5.1 Subcritical flow in a short domain 

In the case of a subcritical flow in a short domain, as in the three that follow, the 
cross section of the channel is rectangular, the bottom is given by the function 
Bi and the length L = 200 m. 
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In this current case, the flow is fixed at inflow and the water height is pre- 
scribed at outflow. We have the following boundary conditions: 

J upstream: g = 20 m'^s"^, 
1 downstream: h = hout- 

The channel is initially dry, with a little puddle downstream, i.e. initial 
conditions are: 

h{x, y) = max(0, hout + ^(200, y) — z{x, y)) and q ~ m'^/s. 
If we take the mean water height 

he^{x) = 0.9 + 0.3cxp 1^-20 

the flow stays subcritical in the whole domain of length L = 200 m. 

3.5.2 Supercritical flow in a short domain 

In this case, the flow and the water height arc fixed at inflow. We have the 
following boundary conditions: 

j upstream: g = 20 m'^s^^ and h = hin, 
\ downstream: free. 

The channel is initially dry, i.e. initial conditions are: 

/i = m, and q = m'^/s. 

If we consider the mean water height 

h^,^{x) = 0.5 + 0.5exp ^-20 

in a channel of length L = 200 m with the Bi shape and vertical boundary, the 
flow is supercritical. 

3.5.3 Smooth transition in a short domain 

In the case of smooth transition in a short domain, the flow is fixed at the inflow. 
We have the following boundary conditions: 

J upstream: g = 20 m'^s"^, 
1 downstream: free. 

The channel is initially dry, i.e. initial conditions are: 

h = m, and q ~ m'^/s. 

The channel is the same as in the previous cases, with a mean water height 
given by 

hex{x) = 1 - O.Stanh ^4 
Under these conditions, the flow is flrst subcritical and becomes supercritical. 
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3.5.4 Hydraulic jump in a short domain 

In the case of a hydraulic jump in a short domain, the flow discharge is fixed 
at the inflow and the water height is prescribed at both inflow and outflow. We 
have the following boundary conditions: 

upstream: g = 20 m'^s"^ and h = hin, 
downstream: h = hout- 

The channel is initially dry, with a little puddle downstream, i.e. initial 
conditions are: 

h{x, y) = max(0, hout + ^(200, y) — z{x, y)) and q = m'^/s. 
We choose the following expression for the mean water height: 

hex{x) = 0.7 + 0.3 (cxp (^^) - l) for m < X < 120 m, 

and 

/ x^x* V 

hex{x) = exp (-p(.T - X*)) ^ fc, ( — J +0(.t) for 120 m < x < 200 m, 

with: 

• X* ^ 120 m, • X** = 200 m, • M = 2, 

• fco = -0.154375, • fci = -0.108189, • fc2 = -2.014310, 
. p = 0.1, . (b{x) = 1.5 exp (o.l - 1 

We obtain a supercritical flow that turns into a subcritical flow through a hy- 
draulic jump. 



3.5.5 Subcritical flow in a long domain 

From now on, the length of the domain is L = 400 m, the boundaries of the 
channel are given by B2 and the cross sections are isoscele trapezoids. 

In this case, the flow is fixed at the inflow and the water height is prescribed 
at the outflow. We have the following boundary conditions: 

J upstream: g = 20 m'^s"^, 
1 downstream: h = hout- 

The channel is initially dry, with a little puddle downstream, i.e. initial 
conditions are: 

h{x, y) = max(0, hout + z(400, y) — z(x, y)) and q = Q w? /&. 
Considering the mean water height 

hM = 0.9 + 0.3 exp (-40 - 1) + 0.2 exp (-35 - , 

the flow is subcritical along the whole channel. 
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3.5.6 Smooth transition followed by a hydraulic jump in a long do- 
main 



In this case, the flow is fixed at the inflow and the water height is prescribed at 
the outflow. We have the foUowing boundary conditions: 

upstream: g = 20 m'^s"^, 
downstream: h = hout- 

The channel is initially dry, with a little puddle downstream, i.e. initial 
conditions are: 

h{x, y) = max(0, hout + ^(400, y) — z{x, y)) and q = m'^/s. 
With the second channel, we define the mean water height by 

= 0.9 + 0.25 (exp (-^) - l) + 0.25exp (l5 - A) 



for m < X < 120 m. 



and 

M 



i=0 



hex(x) ~ exp {-~p{x — X*)) ^ ki ( J +(j){x) for 120 m < a; < 200 m. 



with: 

• x* = 120 m, • X** = 400 m, • Af = 2, 

• fco = -0.183691, • ki = 1.519577, • fca -18.234429, 

. p = 0.09, . cp{x) = 1.5 exp (o.l6 - l) ) - 0.3 exp (2(^-1 

Starting with a subcritical flow, we get a smooth transition to a supercriti- 
cal zone, and, through a hydraulic jump, the flow becomes subcritical in the 
remaining of the domain. 



4 Transitory solutions 

In section 3, we gave steady-state solutions of increasing difficulties. These 
solutions can be used to check if the numerical methods are able to kecp/catch 
steady-state flows. But even if the initial condition differs from the expected 
steady state, we do not have information about the transitory behavior. Thus, 
in this section, we list transitory solutions that may improve the validation of the 
numerical methods. Moreover, as most of these cases have wet/dry transitions, 
one can check the ability of the schemes to capture the evolution of these fronts 
{e.g. some methods may fail and give negative water height). At last, we give 
some periodic transitory solutions in order to check whether the schemes are 
numerically diffusive or not. 

Table 3 lists all transitory solutions available in SWASHES and outlines their 
main features. 
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Transitory solutions 


Mlopc 


l-'riction 




Typo 


Description 


§ 


Reference 


Null 


Const. 


Var- 


Man- 


D.-W. 


Other 


Null 


Comments 


Dam breaks 


Dam break 
on wet domain without friction 


4.1.1 


[46] 


X 












X 


Moving shock. ID 


Dam break 
on dry domain without friction 


4.1.2 


[42] 


X 












X 


Wet-dry transition. ID 


Dam break 
on dry domain with friction 


4.1.3 


[20] 


X 








X 






Wet-dry transition. ID 


Oscillations 
without damping 


Planar surface 
in a parabola 


4.2.1 


[48] 






X 








X 


Wet-dry transition. ID 


Kadialiy-symmctrical 
paraboloid 


4.2,2 


[48] 






X 








X 


Wet-dry transition. 2D 


Planar surface 
in a paraboloid 


4.2,2 


[48] 






X 








X 


Wet-dry transition. 2D 


Oscillations 
with damping 


Planar surface 
in a parabola with friction 


4.2.3 


[43] 






X 






X 




Wet-dry transition. ID 



Canst.: Constant; Var.: Variable; Man.: Manning; D.-W.: Darcy-Weisbach 



Table 3: Analytic solutions for shallow flow equations and their main features 
— Transitory cases 

4.1 Dam breaks 

In this section, we are interested in dam break solutions of increasing complexity 
on a flat topography namely Stoker's, Ritter's and Dressler's solutions. The 
analysis of dam break flow is part of dam design and safety analysis: dam 
breaks can release an enormous amount of water in a short time. This could be 
a threat to human life and to the infrastructures. To quantify the associated 
risk, a detailed description of the dam break flood wave is required. Research on 
dam break started more than a century ago. In 1892, Ritter was the first who 
studied the problem, deriving an analytic solution based on the characteristics 
method (all the following solutions are generalizations of his method) . He gave 
the solution for a dam break on a dry bed without friction (in particular, he 
considered an ideal fluid flow at the wavefront): it gives a parabolic water height 
profile connecting the upstream undisturbed region to the wet /dry transition 
point. In the 1950's, Dressier (see [20]) and Whitham (see [54]) derived analytic 
expressions for dam breaks on a dry bed including the effect of bed resistance 
with Chezy friction law. They both proved that the solution is equal to Ritter's 
solution behind the wave tip. But Dressier neglects the tip region (so his solution 
gives the location of the tip but not its shape) whereas Whitham's approach, 
by treating the tip thanks to an integral method, is more complete. Dressier 
compared these two solutions on experimental data [21]. A few years later. 
Stoker generalized Ritter's solution for a wet bed downstream the dam to avoid 
wet / dry transition. 

Let us mention some other dam break solutions but we do not detail their 
expressions here. Ritter's solution has been generalized to a trapezoidal cross 
section channel thanks to Taylor's series in [55]. Dam break flows are also exam- 
ined for problems in hydraulic or coastal engineering for example to discuss the 
behavior of a strong bore, caused by a tsunami, over a uniformly sloping beach. 
Thus in 1985 Matsutomi gave a solution of a dam break on a uniformly sloping 
bottom (as mentioned in [39]). Another contribution is the one of Chanson, 
who generalized Dressler's and Whitham's dam break solutions to turbulent 
and laminar flows with horizontal and sloping bottom [12, 13]. 

4.1.1 Dam break on a wet domain without friction 

In the shallow water community, Stoker's solution or dam break on a wet 
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domain is a classical case (introduced first in [46, p. 333-341]). This is a classical 
Riemann problem: its analog in compressible gas dynamics is the Sod tube [45] 
and in blood flow dynamics the ideal tourniquet [18]. In this section, we con- 
sider an ideal dam break on a wet domain, i.e. the dam break is instanteneous, 
the bottom is flat and there is no friction. We obtain an analytic solution of 
this case thanks to the characteristics method. 



The initial condition for this configuration is the following Riemann problem 
h{x) = l' 



for m < a: < Xo, 
for Xq < X < L, 



with hi > hr and u[x) = m/s. 



At time t> 0, we have a left- going rarefaction wave (or a part of parabola 
between XA{i) and XB{t)) that reduces the initial depth hi into h„i, and a right- 
going shock (located in xc{t)) that increases the intial height hr into hm- For 
each time t > 0, the analytic solution is given by 



h{t, x) 



where 



4 

J2 



X \ ' 



'ghl - — 



2t) 



u(t, x) 



ni/s 
2 



if a; < XAit), 

2 (^Y + if ^A[i) <x< XB{t), 

2 {Vghi - Cm) if xsit) < X < xc{t), 
m/s if xc{t) < X, 



XA{t) = XQ-t^/ghi, XB{t) = xo+t {2-/ghi - 3c„j^ and xc{t) = xo+t — "'^^2^^^^ — ~ 



withcm = \/ghm solution of -SghrCm"^ {ghi - c™^) -t-(cm^ - gKY 
0. 



Cm - ghr 

'^+ghr) = 



This solution tests whether the code gives the location of the moving shock 
properly. 

In SWASHES, we take the following parameters for the dam: hi = 0.005 m, 
hr = 0.001 m, = 5 m, L = 10 m and T = 6 s. 



4.1.2 Dam break on a dry domain without friction 

Let us now look at Rittcr's solution [42]: this is an ideal dam break (with 
a reservoir of constant height hi) on a dry domain, i.e. as for the Stoker's 
solution, the dam break is instantaneous, the bottom is fiat and there is no 
friction. The initial condition (Riemann problem) is modified and reads: 

, , s _ J /i; > for m < a; < xo, 

1 /ir = m for Xq < X < L, 

and u{x) = m/s. At time t > 0, the free surface is the constant water height 
(hi) at rest connected to a dry zone (hr) by a parabola. This parabola is limited 
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upstream (resp. downstream) by the abscissa XA{t) (resp. xsit)). The analytic 
solution is given by 

if a: < XA{t), 
if XAit) < X < xsit), 
if xsit) < X, 

where 

XA{t) = a;o — t\/^hi and xsit) = xq + 2t\/ ghi, 

This solution shows if the scheme is able to locate and treat correctly the 
wet/dry transition. It also emphasizes whether the scheme preserves the posi- 
tivity of the water height, as this property is usually violated near the wetting 
front. 

In SWASHES, we consider the numerical values: hi = 0.005 m, = 5 m, L = 
10 m and T = 6 s. 




4.1.3 Dressler's dam break with friction 

In this section, we consider a dam break on a dry domain with a friction term 
[20]. In the literature we may find several approaches for this case. Although 
it is not complete in the wave tip (behind the wet-dry transition), we present 
here Dressler's approach. Dressier considered Chczy friction law and used a 
perturbation method in Ritter's method, i.e. u and h are expanded as power 
series in the friction coefficient C/ = 
The initial condition is 



h{x) 



hi >0 
hr = m 



for m < a; < xq, 
for xq < X < L, 



and u{x) = m/s. Dressler's first order developments for the flow resistance 
give the following corrected water height and velocity 



where 



and 



a2 



hco{x.j i) 

9 



[x,t) 



2 _^ .T .g2 



3 

ghi 



2x 



6 



ai = 



12 



5 2- 



2- 



ty/ghi 



ty/ghi 

i 8^3 
! ^ T89 



4^/3 
l35 



ty/ghl 



ty/ghi 



3/2 



7 2 



3/2 



108 



ty/ghi 



(16) 



With this approach, four regions are considered: from upstream to downstream, 
a steady state region ((ft,;,0) for x < x^(i)), a corrected region {{hco,Uco) for 
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XA(t) < X < xrit)), the tip region (for XT{t) < x < xsit)) and the dry region 
((0, 0) for xsit) < x). In the tip region, friction term is preponderant thus (16) 
is no more vahd. In the corrected region, the velocity increases with x. Dressier 
assumed that at xt (t) the velocity reaches the maximum of Uco and that the ve- 
locity is constant in space in the tip region utip{t) = Tii&yixi^^^^(^t),XB{t)] Uco{x,t). 
The analytic solution is then given by 



h{t, x) 



and with 



hi 

m 

XAit) = Xo 



u{t, x) 



m/s 

Uco{x,t) 
Utipit) 

m/s 



and xsit) 



if a; < XA{t), 

if XA{t) < X < XT{t), 

if XT{t) < X < XB(t), 

if XB{t) < X, 



We should remark that with this approach the water height is not modified in 
the tip zone. This is a limit of Dressler's approach. Thus we coded the second 
order interpolation used in [50, 51] (not detailed here) and recommanded by 
Valerio Calcffi^. 

Even if we have no information concerning the shape of the wave tip, this case 
shows if the scheme is able to locate and treat correctly the wet/dry transition. 

In SWASHES, we have /i/ = 6 m, = 1000 m, C = 40 m^/^/s (Chczy 
coefficient), L = 2000 ni and T = 40 s. 



4.2 Oscillations 

In this section, we arc interested in Thacker's and Sampson's solutions. These 
are analytic solutions with a variable slope (in space) for which the wet/dry 
transitions are moving. Such moving-boundary solutions are of great interest 
in communities interested in tsunami run-up and ocean flow simulations (see 
among others [37], [33] and [41]). A prime motivation for these solutions is 
to provide tests for numerical techniques and codes in wet/dry transitions on 
varying topographies. The first moving-boundary solutions of Shallow- Water 
equations for a water wave climbing a linearly sloping beach is obtained in [11] 
(thanks to a hodograph transformation) . Using Shallow- Water equations in La- 
grangian form. Miles and Ball [40] and [3] mentioned exact moving-boundary so- 
lutions in a parabolic trough and in a paraboloid of revolution. In [48] , Thackcr 
shows exact moving boundary solutions using Eulerian equations. His approach 
was first to make assumptions about the nature of the motion and then to solve 
the basin shape in which that motion is possible. His solutions are analytic 
periodic solutions (there is no damping) with Coriolis effect. Some of these so- 
lutions were generalised by Sampson et al. [43, 44] by adding damping due to 
a linear friction term. 

Thacker's solutions described here do not take into account Coriolis effect. In 
the ID case, the topography is a parabola and in the 2D case it is a paraboloid. 
These solutions test the ability of schemes to simulate flows with comings and 

^ Personal communication 
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goings and, as the water height is periodic in time, the numerical diffusion of 
the scheme. Similarly to the dam breaks on a dry domain, Thacker's solutions 
also test wet / dry transition. 



4.2.1 Planar surface in a parabola without friction 

Each solution written by Thacker has two dimensions in space [48]. The exact 
solution described here is a simplification to ID of an artificially 2D Thacker's 
solution. Indeed, for this solution, Thacker considered an infinite channel with 
a parabolic cross section but the velocity has only one nonzero component (or- 
thogonal to the axis of the channel). This case provides us with a relevant test 
in ID for shallow water model because it deals with a sloping bed as well as 
with wetting and drying. The topography is a parabolic bowl given by 



1 



L 



- 1 



and the initial condition on the water height is 

2 



1 



L 



h{x) 



-ho - x~- 



B 



- 1 J for a;i(0) < X < .T2(0), 
m otherwise, 

with B = \/2gho/{2a) and for the velocity u{x) = m/s. Thacker's solution 
is a periodic solution (without friction) and the free surface remains planar in 
time. The analytic solution is 



h{t, x) 




m 



B 



( V^gho 



1 



for xi (t) < X < X2 (t) , 
otherwise, 



u{t, x) 



B sin 



[ a 



m/s 



for xi{t) < X < X2(t), 



otherwise. 



where xi(t) and X2(t) are the locations of wet /dry interfaces at time t 



xi{t) = - -cos 



X2 (i) = - 2 cos 



y/2gho 



In SWASHES, we consider a = 1 ni, ho = 0.5 m, L = 4 m and T = 10.0303 s 
(5 periods). 
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4.2.2 Two dimensional cases 



Several two dimensional exact solutions with moving boundaries were developed 
by Thacker. Most of them include the Coriolis force that we do not consider 
here (for further information, see [48]). These solutions are periodic in time 
with moving wet/dry transitions. They provide perfect tests for shallow water 
as they deal with bed slope and wetting/drying with two dimensional effects. 
Moreover, as the solution is exact without discontinuity, it is very appropriate 
to verify the accuracy of a numerical method. 



Radially-symmetrical paraboloid The two dimensional case presented here 
is a radially symmetrical oscillating paraboloid [48]. The solution is periodic 
without damping (i.e. no friction). The topography is a paraboloid of revolu- 
tion defined by 



z{r)^-ho[l-^] (17) 



with r = \/{x — L/2y + {y — L/2Y for each {x, y) in [0; L] x [0; i], where is 
the water depth at the central point of the domain for a zero elevation and a is 
the distance from this central point to the zero elevation of the shoreline. The 
solution is given by: 




^/l-A^ 
1 — Acos(ajt) 
1 fl 

— v4cos(cj<) 

1 

- ^cos(wt) 




-1 \-z{r) 



where the frequency a; is defined as w = y/Sgho/a, vq is the distance from 
the central point to the point where the shoreline is initially located and A = 
(a^ - r^)/(a2 + r^) (Figure 4). 

The analytic solution at t = s is taken as initial condition. 

In SWASHES, we consider a = 1 m, tq = 0.8 m, ft-o = 0.1 m, L = 4 m and 
T = 3^. 

Planar surface in a paraboloid For this second Thacker's 2D case, the 
moving shoreline is a circle and the topography is again given by (17). The free 
surface has a periodic motion and remains planar in time [48] . To visualize this 
case, one can think of a glass with some liquid in rotation inside. 
The exact periodic solution is given by: 

h{x, y,t) = ^ ^2 (^x - -^^ cos(a;t) + 2 (^J/ - sm{ujt) - - z{x, y) 

u{x,y,t) ~ — ryw sin(ajt) 
v{x, y, t) = rjuj cos(wt) 



for each (x, y) in [0; L] x [0; L], where the frequency w is defined as oj = ^Jigho/a 
and ?7 is a parameter. 

Here again, the analytic solution at < = s is taken as initial condition. 

In SWASHES, we consider a = 1 m, /iq = 0.1 m, r/ = 0.5, i = 4 m and 
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Figure 4: Notations for Thacker's axisymmetrical solution 



4.2.3 Planar surface in a parabola with friction 



Considering a linear friction term {i.e. Sf = Tu/g) in system (2) (with R = 
m/s and ^ = m^/s) with Thacker's approach, Sampson et al. got moving 
boundaries solutions with damping [43, 44]. These solutions provide a set of ID 
benchmarks for numerical techniques in wet/dry transitions on varying topogra- 
phies (as Thacker's solutions) and with a friction term. One of these solutions 
is presented here. The topography is a parabolic bowl given by 

L 



z{x) = ho 

where ho and a are two parameters and x Cz [0, L]. The initial free surface is 
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for a;i(0) < a; < X2(0), 

else, 



with B a constant, s = y^p^ — /2 and p 
remains planar in time 

a^B^e--' 



^/Sghoja^ . The free surface 



{z+h){t,x) 



ho 



Sg^ho 



ST sin(2st) 



s2 cos(2st) 



4.9 



tB 

Bs cos{st) H — — sin(st) 
for Xi{t) < X < X2{t) , 

else 



z{x) 

and the velocity is given by 

Be-''*/'^ sm{st) for xi{t) <x< X2{t), 
Om/s else. 



u{t, x) 
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The wet/dry transitions are located at xi(t) and X2{t) 



xi{t) 



2gho 



( 



Bs cos(.st) 
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+ 
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Bs cos(si) 
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+ 
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2 



In SWASHES, we consider 
5 m/s, L = 10, 000 m and T = 



a = 3,000 ni, ho = 10 m, r = 0.001 s"\ B 
6,000 s. 



5 The SWASHES 



software 



In this section, we describe the Shallow Water Analytic Solutions for Hydraulic 
and Environmental Studies (SWASHES) software. At the moment, SWASHES 
includes all the analytic solutions given in this paper. The source code is 
freely available to the community through the SWASHES repository hosted at 
http://www.univ-orleans.fr/mapmo/soft/SWASHES. It is distributed under 
CeCILL-V2 (GPL-compatible) free software license. 

When running the program, the user must specify in the command line the 
choice of the solution (namely the dimension, the type, the domain and the 
number of the solution) as well as the number of cells for the discretization 
of the analytic solution. The solution is computed and can be redirected in a 
gnuplot-compatible ASCII file. 

SWASHES is written in object-oriented ISO C++. The program is struc- 
tured as follows: each type of solutions, such as bump, dam-break, Thacker, 
etc., is written in a specific class. This structure gives the opportunity to easily 
implement a new solution, whether in a class that already exists (for example 
a new Mac Donald type solution), or in a new class. Each analytic solution is 
coded with specific parameters (most of them taken from [16]). In fact, all the 
parameters are written in the code, except the number of cells. 

We claim that such a library can be useful for developers of Shallow- Water 
codes to evaluate the performances and properties of their own code, each ana- 
lytic solution being a potential piece of benchmark. Depending on the targeted 
applications and considering the wide range of fiow conditions available among 
the analytic solutions, developers may select a subset of the analytic solutions 
available in SWASHES. We recommend not to change the values of the pa- 
rameters to ease the comparison of numerical methods among different codes. 
However, it may be legitimate to modify these parameters to adapt the case to 
other specific requirements (friction coefficient, dam break height, rain intensity, 
etc.). This can be easily done in SWASHES but, in that case, the code must be 
renamed to avoid confusions. 

6 Some numerical results: comparison with FullSWOF 
approximate solutions 

SWASHES was created because we have been developing a software for the 
resolution of Shallow- Water equations, namely FullSWOF, and we wanted to 
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validate it against analytic solutions. To illustrate the use of SWASHES in 
a practical case, we first give a short description of the ID and 2D codes of 
FuUSWOF and then we compare the results of FullSWOF with the analytic so- 
lutions. The comparisons between FullSWOF results and the analytic solutions 
is based on the relative error in percentage of the water height, using, of course, 
the analytic solution as a reference. This percentage is positive when FullSWOF 
overestimates the water height and negative when it underestimates it. 

6.1 The FullSWOF program 

FullSWOF (Full Shallow- Water equations for Overland Flow) is an object- 
oriented C++ code (free software and GPL-compatible license CeCILL-V2^. 
Source code available at http : //www .univ-orleans . f r/mapmo/ sof t/FullSWOF/) 
developed in the framework of the multidisciplinary project METHODE (see 
[15, 16] and http://www.univ-orleans.fr/mapmo/methode/). We briefly de- 
scribe here the principles of the numerical methods used in FullSW0F_2D. The 
main strategy consists in a finite volume method on a structured mesh in two 
space dimensions. Structured meshes have been chosen because on the one 
hand digital topographic maps are often provided on such grids, and, on the 
other hand, it allows to develop numerical schemes in one space dimension 
(implemented in FullSWOF_lD), the extension to dimension two being then 
straightforward. Finite volume method ensures by construction the conserva- 
tion of the water mass, and is coupled with the hydrostatic reconstruction [2, 7] 
to deal with the topography source term. This reconstruction preserves the 
positivity of the water height and provides a well-balanced scheme (notion in- 
troduced in [28]) i.e. it preserves at least hydrostatic equilibrium (9) (typically 
puddles and lakes). Several numerical fluxes and second order reconstructions 
are implemented. Currently, we recommend, based on [16], to use the second 
order scheme with MUSCL reconstruction [52] and HLL flux [29]. FuhSWOF is 
structured in order to ease the implementation of new numerical methods. 

6.2 Examples in one dimension 

In this part, wc give the results obtained with FullSWOF_lD for three one- 
dimensional cases. For these examples, we have 500 cells on the FullSWOF_lD 
domain but 2500 for semi- analytic solutions (see Remark 1). In the flrst two 
flgurcs, wc also plotted the critical height, in order to show directly whether the 
flow is fluvial or torrential. 

6.2.1 Transcritical flow with shock 

On Figure 5, we plotted the solution for a transcritical flow with shock (Sec- 
tion 3.1.5) for a time large enough to attain the steady state, namely t = 100 s. 

Overall, the numerical result of FullSWOF_lD and the analytic solution are 
in very good agreement. By looking carefully at the differences in water height, 
it appears that the difference is extremely low before the bump (—4 x 10^^ m, 
i.e. -0.001%). Right at the beginning of the bump, this difference increases 
to +0.06%, but on a single cell. Over the bump, differences are alternatively 

^http: / /www. cecill. info/index. en. html 
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Figure 5: Results of FullSWOF_lD for a transcritical flow with shock 



positive and negative which leads, overall, to a good estimate. A maximum 
(+1.2%) is reached exactly at the top of the bump {i.e. at the transition from 
subcritical flow to supercritical flow). Just before the shock, the difference is 
only +0.03%. The largest difference in height (+100%) is achieved at the shock 
and affects a single cell. While for the analytic solution the shock is extremely 
sharp, in FullSWOF_lD it spans over four cells. After the shock and up to the 
outlet, the heights computed by FullSWOF_lD remain lower than the heights of 
the analytic solution, going from -1% (after the shock) to -0.01% (at the outflow 
boundary condition). 

6.2.2 Smooth transition and shock in a short domain 

Figure 6 shows the case of a smooth transition and a shock in a short domain 
(Section 3.2.2). The final time is t = 150 s. 

FullSWOF_lD result is very close to the analytic solution. From a; = ni to 
the subcritical-to-supercritical transition, FullSWOF_lD undcrcstimes slightly 
the water height and the difference grows smoothly with x (to -0.1%). Around 
a; = 45 m (i.e. close to the subcritical-to-supercritical transition), the difference 
remains negative but oscillates, reaching a maximum difference of -0.22%. After 
this transition, FullSWOF_lD continues to underestimate the water height and 
this difference grows smoothly up to -0.5% right before the shock. As for the 
comparison with the transcritical flow with shock (see Section 6.2.1 ), the max- 
imum difference is reached exactly at the shock {x = 66.6 m): FullSWOF_lD 
overestimates water height by +24%. Right after the shock, the overestimation 
is only +1% and decreases continuously downstream, reaching 0% at the outlet. 

6.2.3 Dam break on a dry domain 

The last figure for the one-dimensional case is Figure 7, with the solution of 
a dam break on a dry domain (Section 4.1.2). We chose the final time equal 
to t = 6 s. Here too, the numerical result of FullSWOF_lD matches with the 
analytic solution well. The height differences on the plateau (h ~ 0.005 m) are 
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Figure 6: Results of FullSWOF.lD for the Mac Donald's type solution with 
a smooth transition and a shock in a short domain, with Manning's friction 
coefficient 
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Figure 7: Results of FullSWOF_lD for a dam break on a dry domain without 
friction 
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null. On this domain, the water is not flowing yet. This shows FullSWOF_lD 
preserves steady-states at rest properly. The analytic solution predicts a kink 
at a; = 3.7 m. In FullSWOF_lD result, it does not show a sharp angle but a 
curve: FullSWOF_lD slightly underestimates the water height (by a maximum 
of -2%) between x = 3.4 m and a; = 4 m. Between a; = 4 m and x = 6.9 m, 
FuUSWOF_lD overestimates constantly the water height (by up to +2%). On 
all the domain between a; = 6.9 m and a; = 7.6 m. the water height predicted by 
FullSWOF_lD is lower than the water height computed by the analytic solution 
(by as much as -100%). At this water front, the surface is wet up to a; = 7.6 m 
according to the analytic solution, while FullSWOF_lD predicts a dry surface for 
X > 7.3 m. It shows the water front predicted by FullSWOF_lD moves a bit too 
slowly. This difference can be due to the degeneracy of the system for vanishing 
water height. For x > 7.6 m, both the analytic solution and FullSWOF_lD 
result give a null water height: on dry areas, FullSWOF_lD does not predict 
positive or negative water heights, and the dry-to-wet transition does not show 
any spurious point, contrarily to other codes. 

6.3 Examples in two dimensions 

We now consider the results given by the two-dimensional software FullSW0F_2D. 

6.3.1 Planar surface in a paraboloid 

In the case of Thacker's planar surface in a paraboloid (see Section 4.2.2), 
FullSW0F_2D was run for three periods on a domain of [0 m; 4 m] x [0 m; 4 m] 
with 500 X 500 cells. To analyse the performances of FullSW0F_2D, wc consider 
a cross section along x (Figure 8). 

Overall, FullSW0F_2D produces a good approximation of the analytic so- 
lution: while the maximum water height is 0.1 m, errors are in the domain 
[—1.55 X 10"^ m; +7.65 x 10^^ m]. At the wet-dry transition located at a; = 
1.5 m, FullSW0F_2D overestimates the water height by about -1-6 x 10~* m and, 
on a single cell, FullSW0F_2D predicts water while the analytic solution gives 
a dry surface. Close to a; = 1.5 m, the overestimation of the water height is up 
to -f43%, but decreases quickly and is always less than -1-5% for x > 1.57 m. 
The overestimation persists up to a; = 2.65 m and then becomes an underes- 
timation. The underestimation tends to grow up to the wet-dry transition at 
X = 3.5 m (reaching -5% at a: = 3.38 m). Starting exactly at x ~ 3.5 m (and on 
the same cell) both FullSW0F_2D and the analytic solution predict no water. 
However, on the cell just before, FullSW0F_2D underestimates the water height 
by -97% (but at this point, the water height given by the analytic solution is 
only 1.6 X lO^'^ m). 

6.3.2 Mac Donald pseudo-2D solutions 

We consider two Mac Donald pseudo-2D solutions. Since FullSW0F_2D does 
not solve the pseudo-2D Shallow- Water system but the full Shallow- Water sys- 
tem in 2D (1), more significant differences are expected. In both cases, FullSW0F_2D 
was run long enough to reach steady-state. 
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Figure 8: Cross section along x of FullSW0F_2D result for Thacker's planar 
surface in a paraboloid 



Supercritical flow in a short domain The case of a supercritical flow in 
a short domain (Section 3.5.2) is computed with 400 x 201 cells on the domain 
[0 m; 200 m] x [0 m; 9.589575 m] , with the topography Bi . The y-averaged result 
of FullSW0F_2D differs from the analytic solution mainly around a; = 100 m 
(Figure 9), with an underestimation of the water height of up to -0.018 m (- 
11.85%). This underestimation occurs on the whole domain and gets closer to 
zero near both the upper and lower boundaries. 



Subcritical flow in a long domain The case of a subcritical flow in a 
long domain (Section 3.5.5) is computed with 800 x 201 cells on the domain 
[0 m; 400 m] x [0 m; 9.98067 m], with the topography i?2- Comparison between 
the y-averaged FullSW0F_2D result and the analytic solution shows clear dif- 
ferences (Figure 10), even if the overall shape of the free surface given by 
FullSW0F_2D matches the analytic solution. FullSW0F_2D underestimates 
the water height on most of the domain. This difference can be up to -0.088 m 
(-8.8%) at a; = 66 m. FuUSW0F_2D overestimates water height for x > 297 m 
and this overestimation can reach up to 4-0.04 m (+4.1%) at x = 334 m. 
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